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The effect of carrier multiplication (CM) in semiconductor nanocrystals is systematically treated 
by employing an exciton scattering approach. Using projection operators, we reduce the Coulomb 
coupled multi-exciton dynamics to scattering dynamics in the space spanning both single- and bi- 
exciton states. We derive a closed set of equations determining the scattering matrix elements. This 
allows us to interpret CM dynamics as a series of odd-order interband scattering events. Using the 
time-dependent density matrix formalism, we provide a rigorous description of the CM dynamics 
induced by a finite-time pump pulse. Within this approach, both processes of single- and bi-exciton 
photogeneration and the consequent population relaxation are treated on the same footing. This 
approach provides a framework for numerical calculations and for comparisons of the quantum 
efficiencies associated with each process. For applications, the limit of weak interband Coulomb 
coupling is considered. Finally, we demonstrate that three previously used theoretical models can 
be recovered as limiting cases of our exciton scattering model. 
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I. INTRODUCTION 

Carrier multiplication (CM) in semiconductor mate- 
rials is the process of more than one electron-hole pair 
generation per single absorbed photon. Here, we con- 
sider the general case in which high energy electron-hole 
pairs consist of free carriers (as typically occurs in bulk 
semiconductors). We also consider the case in which the 
carriers are confined exciton states as occurs in semi- 
conductor nanocrystals (NCs). CM is naturally charac- 
terized by the related Quantum Efficiency (QE) which is 
the number of electron-hole pairs generated per absorbed 
photon. CM is also characterized by the activation en- 
ergy threshold (AET) below which CM becomes negligi- 
ble. Extensive studies of CM are motivated by potential 
applications in photovoltaic, photoelectrochemical, and 
energy storage devices .1^^ 

CM was first investigated in bulk materials using pho- 
tocurrent measurement s,'^^^ and was recently revisited 
using terahertz time-domaing spectroscopy.^^ The theory 
of CM in bulk treats CM as a sequence of the primary 
photoexcitation event in which a single electron-hole pair 
is created by a photon, and the secondary process of 
the electron and hole population relaxation during which 
CM occurs .2^^ The population relaxation dynamics is 
a competition between the impact ionization process in 
which the excess kinetic energy of the hot electron or hole 
is transferred to create another electron- ho lepaii^^f^^l^ 
and the process of phonon-assisted cooling .fl^^ 

In bulk, strict energy and quasi-momentumconserva- 
tion constraints determine the values of and QE 

as a function of the absorbed photon eneigy^^. For a va- 
riety of semiconductor materials, the lower boundary of 
AET is fou nd to be about 3Eg where E^is the bulk band 
gap energy.f^Sl^ However, photo currentP^ and opticaP^ 
measurements have demonstrated that the AET for most 



materials is > AEgW^ 

In semiconductor NCs, it is expected that the following 
three processes should lead to an increase in QE and a 
decrease in AET: relaxation of the quasi- momentum con- 
servation constraint!^, a decrease in the phonon-assisted 
relaxation rat^^, and an enhancement of Coulomb inter- 
action between the carrierJ^^^. Efficient CM has been 
reported in colloidal NCs using time-resolved transient 
absorption (pump-probe) and time-resolved photolumi- 
nescence techniques. ^^ ^^ Reported values of AET vary 
in the range of 2 — 3Eg with Eg being the NCs band gap 
energy. By using a bulk-type model with relaxed quasi- 
momentum conservation rule, it has been further specu- 
lated that, depending on the ratio of electron and hole 
effective masses, the AET can reach a minimum value of 
2Eg satisfying the energy conservation constraint. 

These experimental results have been challenged by a 
number of reports claiming significantly lower QE and 
even the absence of the CM effect. CM has been 
further reconfirmed, however, with observed values of 
QE varying in a broad range starting below the QE 
in bulk materials.^^^^' The variation of QE co uld pos- 
sibly arise from experimental inaccuracies sample- 
to-sample variation in surface preparation j ^^ l ^^ l and ex- 
traneous effects such as photochargin^^. These is- 
sues raise an important question: What are the specific 
quantum-confinement-induced features that distinguish 
CM in NCs from CM in bulk semiconductors?^^ Address- 
ing this question requires theoretical insight. Currently, 
there are three separate models outlined below proposing 
different mechanisms for CM in NC: the Coherent Super- 
position Model, the Direct Photogeneration Model, and 
the Impact Ionization Model. We describe each of these 
in the following few paragraphs. 

The Coherent Superposition Model of resonant (almost 
degenerate) single- and bi-exciton states is based on the 
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density matrix formalism. It was proposed by Shabaev, 
Efros, and Nozik.^^ This model states that, in contrast 
to bulk materials, the primary event of single photon 
absorption in NCs leads to the preparation of coherent 
superpositions (oscillations) between the single- and bi- 
exciton states that are almost degenerate. The secondary 
process of phonon-induced intraband relaxation merely 
stabilizes the populations leading to efficient bi-exciton 
production due to the fast bi-exciton intraband relax- 
ation rate. No experimental observations of these os- 
cillations have been reported yet. The enhancement of 
QE according to this model requires a strong Coulomb 
coupling between single- and bi-exciton states. This en- 
hancement has not been confirmed experimentally. This 
model ignores the effects of the single- /bi-exciton density 
of states (DOS) by considering only one single-exciton 
and one bi-exciton states coupled through Coulomb in- 
teractions. As we demonstrate in this paper, this model 
also misses the CM pathway that involves the phonon- 
assisted relaxation channel between single- and bi-exciton 
states. 

Assuming weak Coulomb coupling between single- and 
bi-exciton states and assuming optical pulse duration 
larger than the dephasing t ime, t he QE can be evaluated 
using Fermi's Golden Rule.'^^'^ This approach, referred 
as the Direct Photogeneration Model, predicts two path- 
ways for direct bi-exciton production during the primary 
photon absorption event. The first pathway, introduced 
by Schaller, Agranovich and Klimov, describes resonant 
bi-exciton generation via virtual single-exciton states. 
The second pathway, considered by Rupasov and Klimov, 
accounts for the non- vanishing Coulomb matrix elements 
between the exciton vacuum (filled valence band) and bi- 
exciton states. This coupling leads to the stabilization of 
bi-exciton populations through resonant intraband opti- 
cal transitions."^^ 

These authors estimate the contributions of their re- 
spective pathways and claim that their pathways become 
efficient in NCs because the quasi-momentum conserva- 
tion constraint is relaxed. The actual enhancement of QE 
comes from the increased bi-exciton DOS compared with 
the single-exciton DOS. Independent quantum chemistry 
calculations confirm the possibility of direct carrier pho- 
togeneration in semiconductor clusters. The drawback 
of the Direct Photogeneration Model is that no secondary 
events of population relaxation on QE are considered. 
In our paper, we will also demonstrate that the addi- 
tional channel associated with the direct excitation of 
single-exciton states and their further scattering to the 
bi-exciton manifold during the interaction with the opti- 
cal pulse as well as the interference of all the pathways 
must be included in the weak Coulomb limit. 

A number of reported calculations suggest that, in con- 
trast to the mechanisms outlined above and similar to 
the bulk materials, CM in NCs occurs solely due to the 
competing phonon-assisted relaxation and impact ion- 
ization processes that follow the primary single-exciton 
photoexcitation event. We will refer to this approach 



as the Impact Ionization Model model throughout this 
paper. Specifically, Franceschetti, An, and Zunger have 
considered the spectral dependence of the impact ion- 
ization rate and Auger recombination (the inverse pro- 
cess) rate using atomistic pseudopotential calculations.!^ 
Allan and Delerue used a tight-binding model to simu- 
late the competing processes of impact ionization and 
phonon-assisted relaxation. Their analysis based on 
their models emphasizes the importance of the high ratio 
of bi- to single-exciton DOS for efficient CM. 

Further development of this approach led to a DOS- 
based comparison of QE due to impact ionization and 
direct photogeneration,^ to evaluation of the band- 
structure effects on QE in a variety of NC materials,!^ 
and to modeling the influence of surface defects on QE.^ 
Interestingly, Rabani and Baer emphasized the impor- 
tance of the trion DOS (in contrast to the bi-exciton 
DOS) directly entering the impact ionization and Auger 
recombination rates, where strict selection rules enter 
through the Coulomb matrix elements.'^ 

Currently, the Coherent Superposition Model, the Di- 
rect Photogeneration Model, and the Impact Ionization 
Model are considered as alternative approaches whose 
applicability is still being debated. We propose a more 
general approach capable of treating the CM dynamics 
in both NCs and bulk materials by accounting for both 
the photogeneration event induced by a finite-time op- 
tical pulse and the population relaxation dynamics, on 
the same footing. This approach can be used to interpo- 
late between strong and weak Coulomb coupling regimes. 
This interpolation can be achieved by treating Coulomb 
interactions between the carriers as multiple-scattering 
events. We achieve this more general approach by in- 
tegrating the scattering theory with the density matrix 
formalism, and we call this approach the Exciton Scat- 
tering Model. 

As a validation of our Exciton Scattering Model, we 
demonstrate that the previously proposed models can 
be recovered as limiting cases, and they are fundamen- 
tally related to each other. We also demonstrate that 
the proposed model predicts additional contributions to 
the Coherent Superposition and the Direct Photogener- 
ation models which have not been considered before. We 
use our Exciton Scattering Model to formulate a closed 
computational scheme for determining QE and ATE. The 
results of the numerical simulations using this approach 
applied to specific semiconductor materials will be re- 
ported in a separate paper. 

This paper is organized as follows: In Sec. |TI| the 
projection operator technique is employed to reduce the 
coupled multi-exciton dynamics to single- and bi-exciton 
scattering dynamics in Hilbert space. We use this tech- 
nique to derive a closed set of equations for the interband 
scattering matrix. In Sec. |III[ we use the density matrix 
formalism combined with a modified exciton scattering 
approach to obtain general expressions for the QE that 
naturally describe the primary event of single- and bi- 
exciton photogeneration due to both a finite-time pump 



pulse and population relaxation dynamics. For numeri- 
cal calculations, the limiting case of weak Coulomb cou- 
pling is introduced in Sec. |IV| in which a closed set of 
equations for the limiting QE is presented. In Sec. [V| 
we discuss connections of the Exciton Scattering Model 
with previously proposed models. Finally, we present our 
conclusions in Sec. IVII 



II. MILTI-EXCITON DYNAMICS IN HILBERT 
SPACE 

In this section, we begin our analysis by introduc- 
ing the many-body electronic Hamiltonian in the multi- 
exciton representation accounting for the contributions 
of all of the Coulomb terms. These terms can be parti- 
tioned into those terms that conserve the total number 
of excitons (and determine their binding energies), and 
those terms that do not conserve the number of excitons, 
giving rise to the CM dynamics. An exact treatment of 
the dynamics of total multi-exciton space is not feasi- 
ble. Therefore, we restrict our dynamics to the reduced 
space spanning single- and bi-exciton states by using the 
projection operator technique. This approach allows us 
to include some of the effects of higher- multiplicity (tri-, 
four-, etc.) exciton states in the dynamics in the re- 
duced space. Since the projected dynamics is restricted 
to coupled single- and bi-exciton manifolds only, it can be 
treated by performing an exact summation of the pertur- 
bation series, in which the odd- order inteibdind scattering 
events describe CM dynamics. 
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FIG. 1: The multi-exciton Hamiltonian in the block-matrix 
representation. The (0|7^|0)-block is the exciton vacuum, and 
the rest of the diagonal (red) blocks are the Hamiltonian com- 
ponents in the single-exciton, bi-exciton, etc. subspaces. The 
green off-diagonal blocks describe the Coulomb interactions 
between the latter components changing exciton multiplic- 
ity by one, and the blue blocks by two. 0-blocks indicate 
null-matrices. Four quadrants separated by the black dashs 
describe the partitioning of the multi-exciton Hamiltonian by 
the projection operators P and Q. The upper left quadrant is 
the projected Hamiltonian (Eqs. (|6| - ([s])) acting in the space 
spanning the single- and bi-exciton states. 



A. The multi-exciton Hamiltonian 

Let us consider the valence and conduction bands of 
a semiconductor NC in which single electron and hole 
states are known within the Hartree-Fock approxima- 
tion (or equivalently within the effective mass envelope 
function formalism). Our many-body electronic Hamil- 
tonian, accounts for these non-interacting single par- 
ticle states, and all possible Coulomb interactions among 
them.E^ An explicit form of this Hamiltonian is given in 
Appendix |a] Not all the Coulomb interaction terms in ?i 
conserve the total number of electrons and holes. How- 
ever, this Hamiltonian does conserve total charge. Con- 
sequently, the dynamics of electrically neutral electron- 
hole pairs (excitons) is uncoupled from the dynamics of 
the charged states. This allows us to focus on the dy- 
namics determined by the multi-exciton Hamiltonian, 

HmX = H^MX + ^MX + ^MX' (1) 

whose derivation is provided in Appendix [Aj 

The block-matrix representation of this multi-exciton 
Hamiltonian is shown in Fie; [l] in which the (0|7i|0) -block 
denotes the exciton vacuum, i.e. the filled valence band, 
with its energy set to zero. The remaining diagonal (red) 



blocks describe the single-exciton, bi-exciton, etc.^ sub- 
spaces. These terms are: 

n>l p>l 

in which \x^) denotes the p-th exciton state with multi- 
plicity, n, and energy, hw^. This energy already includes 
the n-particle binding interactions which can be calcu- 
lated, e.g. by block-wise matrix diagonalization. 

The off-diagonal Coulomb interaction (green and blue) 
blocks do not conserve the total number of electrons and 
holes. They describe the interactions between exciton 
states with different multiplicity. Specifically, the green 
off-diagonal blocks describe the processes changing mul- 
tiplicity by one. The blue off-diagonal blocks change 
multiplicity by two. Note that the bi-exciton states are 
coupled to the vacuum, whereas the single-exciton states 
are uncoupled from the vacuum. This is a result of the 
Hartree-Fock representation eliminating the latter inter- 
actions. The general expression for these off-diagonal 
terms in the multi-exciton Hamiltonian is 

n pq 
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in which F^^f +^ = (x^|7i|x^+^) is the interband multi- 
exciton interaction matrix element with z = 1,2 describ- 
ing the multiplicity variation. 



Projected dynamics in single- and bi-exciton 
space 



couplings, V^q'^. Explicit representations for interaction 
matrix elements in terms of the single-particle couplings, 
and related matrix equations defining the single- and bi- 
exciton states are provided in Appendix [A} 

The dynamics restricted to the subspace of interest is 
fully defined by the projected propagator G{t) = PU {t)P 
whose representation in the frequency domain i^^ 



The dynamics in the total multi-exciton Hilbert space 
is fully defined by the propagator: 



U{t) = e(t)exp 



-ih ^HMxt 



(4) 



whose calculation and general analysis is not feasible due 
to the rapidly growing number of multi-exciton states. 
Thus, we introduce a convenient representation allowing 
us to approximate calculations of this propagator. 

The CM processes excited near the AET should pri- 
marily result in photogeneration of single- and bi-exciton 
states. On the other hand, the higher-multiplicity exci- 
ton states can still affect their dynamics. Therefore, we 
consider a dynamics projected onto the space spanned 
by the single- and bi-exciton states, and seek the con- 
ditions allowing us to neglect the effects of the higher- 
multiplicity states. This can be naturally done by intro- 
ducing the following projection operator onto the space 
of single- and bi-exciton states. 



P = ^ \Xa){Xa \ + ^ \xXk){xXk\, 



(5) 



a>0 



k>l 



and the complementary projection operator, Q = I — P 
onto the rest of the multi-exciton space. / denotes the 
identity operator in total multi-exciton space. To distin- 
guish the single- and bi-exciton states from the remaining 
higher-multiplicity states, \Xp) where n > 3, we have in- 
troduced their new notations \xa) and \xxk)^ respectively, 
and use \xo) to denote the exciton vacuum. This notation 
will be used through out this paper. 

The partitioning of the multi-exciton Hamiltonian bv 
the projection operators, P and Q, is illustrated in Fig.[lj 
where the projected Hamiltonian H = PHmxP is the 
sum of two terms 

H = Ho + Vc- (6) 

Here, the first term. 

Ho = '^\Xa)hj^{Xa\ ^'^\xXk)hj^''{xXk\, (7) 



a>l 



k>l 



describes non-interacting single- and bi-exciton states, 
and the second term 

a>0 k>l 

represents the interband Coulomb interactions, V^'^^^ 
between the states, as well as the vacuum to bi-exciton 



G(u;) 



w - heffioj) +i7 



(9) 



where 7 is the finite broadening associated with the 
exciton-phonon coupling. The non-local effective Hamil- 
tonian entering this Green function can be partitioned 
into the sum of the diagonal and off-diagonal terms 



heff{oj) = h{uj) + v{uj), 
which have the following forms 



(10) 



(11) 

(12) 



respectively. The first terms in Eqs. ( pTj ) and ( p!2j ) 
are components of the projected Hamiltonian h~^H 
(Eq. (|6|-([8|) describing the propagation of the coupled 
single- and bi-exciton states. The second terms, account- 
ing for the effect of the higher- multiplicity exciton states, 
are the diagonal, kd{oj)^ and off-diagonal, ko{oj)^ compo- 
nents of the non-local memory kernel, respectively. 

The memory kernel components can be explicitly rep- 
resented in the multi-exciton bases as 



kd{uj) h 



E 



El 

nm>3 pq^:! 



(13) 



ab>0 



ko{uj) 



kl>l 



E EEEi 

nrfi>3pq>l a>0 k>l 



(14) 



where V^^^^ i^i^q'^) ^1^^ interaction matrix elements 
(Eq. ([3|) which couple single-exciton (bi-exciton) states 
with the states of multiplicity nfh > 3. The matrix ele- 
ments G^'^(co') = {Xp\G{uj)\x^) of the propagator 



G{oj) 



CO - h ^QHmxQ 



(15) 



describe the projected dynamics in the higher- 
multiplicity exciton space defined by Q (the lower right 
quadrant in Fig[T]). 
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FIG. 2: Feynman diagram representation of the projected 
propagator, (5, in terms of the scattering operator, T. Panels 
(a) and (b) show the single- and bi-exciton propagator com- 
ponents, respectively, (c) The interband component of the 
propagator mixing the single- and bi-exciton states. 



(a) ia^ ia^ g"" ia^ 

= + 00 + 



(b) J^xx j^^xx ^^xx^xx^^xx ^^xx^xx^^xx ^xx^^xx 

(c) fx,xx _-^X,XX _-fjX,XX^XX ^^xx _-^x,xx^xx -XX ^gXX ^^xx 

<2- i , Oh + O^hOh i.. 

FIG. 3: Feynman diagram expansion for the scattering ma- 
trix, T. Panels (a) and (b) describe even-order scatter- 
ing events contributing to the single-exciton and bi-exciton 
scattering matrix respectively, (c) The odd-order scattering 
events changing the exciton multiplicity contribute to the in- 
terband scattering matrix. The latter processes give rise to 
CM. 



The use of projection operators allows us to map the 
propagator acting in the multi-exciton space (Eq. Q) to 
the projected propagator acting in the space of single- 
and bi-exciton states (Eqs. (|9|-([T5|). This representa- 
tion is exact, since the effect of the higher-multiplicity 
exciton states is fully accounted for through the memory 
kernel (Eqs. (13) and ([l4|). The dynamics of interest can 
now be interpreted as the uncoupled propagation within 
single- and bi-exciton manifolds described by the zeroth- 
order Green function, 



g{uj) = i uj — h{(jj) + Z7 



(16) 



and the scattering events between these manifolds in- 
duced by the interaction operator v{uo). 



C. Single- and bi-exciton scattering model 

To apply the scattering matrix formalism, we represent 
the projected Green function (Eq. (|9|) as a 2 x 2 block 
matrix 



The Fourier transformation of Eq. (pTl) 



(17) 



G{t)= r ^G{uj)exp{-icot), (18) 
J-co 27r 

defines time-evolution of tfie single- and bi-exciton states 

\xa{t)) = Y,G:,{t)\xm)+Y.^::Tit)\xxm),{m 

b>l k>l 

\xx,{t)) = J2Gl^f{t)K{0)) + ^Glf{t)\xx^m20) 



a>0 



1>1 



According to Eqs. ([19| and ( |2Q| ), the matrix elements 
G^^it) and G^f{t) associated with the diagonal blocks in 



Eq. (17) determine the intraband propagation, and the 



matrix elements G^'^^{t) associated with the off-diagonal 



blocks in Eq. (17) describe the interband scattering pro- 
cesses mixing the single- and bi-exciton states. 

Within the scattering matrix formalism, the propaga- 
tor, G, satisfies the following equation?^ 



G{uj) = g{uj) + g{uj)f{uj)g{uj) 



(21) 



in which g{uj) is the intraband zeroth-order Green func- 
tion introduced in Eq. (16). In the modified block-matrix 



representation, this Green function is 

»(")=(''T'»4))- 

Here, the diagonal blocks, g^{uj) and g^^{uo)^ can be de- 
termined numerically using Eq. (16) with matrix inver- 
sion. Finally, the scattering operator in the same repre- 
sentation is 



T 



^(w) f^^ 



(CO) 



(23) 



containing both single-exciton (bi-exciton) component, 
rpx f^rpxx^^ interband components, T^'^^. The Feyn- 
man diagram representation of Eq. (21 ) is given in Fig.[2j 
To find the solution of Eqs. (21)-(23), we need to know 



Eq. (23). 



the form of the matrix elements of 

To obtain a closed set of equations for the scattering 
matrix T, the projected propagator (Eq. ^) should be 
expanded in a power series of the interband coupling op- 
erator, v{uj) (Eq. (12)). These expansion terms can be 



further regrouped to match the form of Eq. (21), leading 



to the diagrammatic expansion of the scattering operator 
shown in Fig. [3] 

In Fig. [sj panels (a) and (b), represent the diagrams 
contributing to the diagonal scattering matrix blocks 
T^{uj) and T^^(a;), respectively. Each term there con- 
tains an even number of vertices, reflecting the even num- 
ber of interband scattering events. This leads to conser- 
vation of the excitons multiplicity, and to renormalization 
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of their energies. The summation of this diagrammatic 
series in panels (a) and (b) results in a set of linear equa- 
tions for the single-exciton (n = x) and the bi-exciton 
(n = xx) scattering matrix elements: 

E [Sif^^ij - Kk{oj)gM] Tl^ioj) = i<T%{o^). (24) 

kl 

Here, the self-energy matrix elements renormalize the 
bare single- and bi-exciton energies, and according to 
Fig. [3] (a) and (b) can be represented as 



n / \ -fc — 2 \ ^ n,m/ \ m/ \ m.n/ \ 

aij{uj)=ih 2^v.'^ {oj)gki{uj)viJ (w), 



(25) 



kl 



where m = xx {m = x) if n = x {n = xx). 

According to Fig. [31(c), the interband scattering ma- 
trix depends on T^^ and can be calculated from the fol- 
lowing linear transformation: 



Ta^l (W) 



mn 



-<'^?r('^)5™„(^)T„T(u;) 



(26) 



This scattering matrix accounts for the odd-oidei scatter- 
ing events changing the multiplicity of the initial exciton 
state, and, therefore, describes the CM dynamics. 

Equations (21)-(26), are exact, since they account 



for all terms entering the multi-exciton Hamiltonian 
(Eqs. ([T])-(3])). Specifically, these terms determine the 
single- and bi-exciton binding energies and the inter- 
band interactions including the effects of the higher- 
multiplicity exciton states. In terms of the diagrammatic 
expansions shown in Figs. [2] and [3) the single and dou- 
ble lines associated with the components of g{uj) and the 
vertices v{uj) are dressed by these interactions. In prac- 
tice, however, only approximate representations for g{uj) 
and v{uj) could be found. For instance, the multi-exciton 
binding energies can be determined approximately or 
even neglected. 

The more difficult task is the evaluation of the mem- 
ory kernel entering the effective Hamiltonian (Eq. ([lO|), 
since the propagator G{uj) (Eq. (15)) cannot be calcu- 
lated exactly. However, the kernel can be calculated ap- 
proximately if the tri-exciton states are accounted for 
only. In this case the kernel will renormalize the single- 
and bi-exucotin resonances showing their hybridization 
with the tri-exciton ones. If the tri-exciton (and higher- 
multiplicity exciton) poles are well-separated from the 
single- and bi-exciton resonances participating in the 
photoexcited dynamics, then the memory kernel can be 
dropped from the projected propagator, G{(jj) (Eq. 
This situation is expected to take place in the vicinity 
of the AET, depending on the strength of the Coulomb 
couplings, V^^^^ and V^^^'^. 

At this point, we focus on the photoinduced dynam- 
ics in the vicinity of the AET only, and for the rest of 
the paper we assume that the memory kernel effect is 
negligible. Therefore, the scattering operator and the 




FIG. 4: Level diagram of CM dynamics in ensemble of NCs. 
(a) Photoexcitation by a pump pulse with central frequency, 
(jjpm, and finite spectral widths results in the generation of 
single- and bi-exciton populations in which the Coulomb scat- 
tering mixes all int erband and intraband dipole transitions 
present in Eq. (28). (b) During the population relaxation, 



both the intraband and the interband processes are mixture of 
the phonon- assisted processes and Coulomb scattering events. 



projected propagator can now be calculated by solving 
the set of linear Eqs. (l24|)-(26) where the zeroth-order 



propagator g{uj) and interband interaction operator v{(jj) 
depend on the projected Hamiltonian (Eqs. ([6|-([8|) only 
and Eqs. ([9]), ^jl, and ([iS])), respectively. 



III. 



PHOTOINDUCED DYNAMICS IN 
LIOUVILLE SPACE 



In this section, we consider the carrier dynamics in an 
ensemble of NCs excited by a pump pulse whose fluence is 
adjusted so that no more than a single photon is absorbed 
per NC. This results in the preparation of no more than 
one single- or bi-exciton state in each NC interacting with 
photons leading to a total population produced by the 
pulse which can be determined by the ensemble average. 

The photoinduced ensemble dynamics is illustrated in 
Fig. [4j Panel (a) shows the exciton photogeneration 
event which occurs on the pump timescale ranging be- 
tween 50 — 100 fs. During the photogeneration, the rel- 
ative number of single- and bi-exciton states produced 
by the pump is determined by the interband scatter- 
ing processes. The photogenerated populations further 
relax on the timescale of 1 — 10 ps as shown in panel 
(b). As we demonstrate below, this relaxation includes 
phonon-assisted cooling to the bottom of the single- and 
bi-exciton bands mixed with the interband population 
transfer due to the Coulomb scattering. The population 
from the bottom of the bi-exciton band finally decays to 
the lowest single-exciton states through Auger recombi- 
nation within > 10 ps. This process (not shown in Fig.[4| 
is typically employed for the experimental determination 
of the bi-exciton production yield, and has no contribu- 
tion to QE. Therefore, we do not consider this process in 
this paper. 
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To include the interaction with the optical pump, we 
extend the projected Hamiltonian as 



Hopt = H + V{t) 
where the following time-dependent term is added 



(27) 



V{t) 



-m 



2b>0 



kl>l 



2k>l 



Here, the optical pulse 



£pra{t) exp{-iujt) 



(29) 



is characterized by the absolute value of the envelope 
function, Spm{t)^ with the widths, r^m, describing the 
pulse duration, and the central frequency. 
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The 

pump envelope and spatial phases do not contribute to 
the population dynamics and therefore are dropped. De- 
tails of the derivation of Hopt are given in Appendix [c) 



According to Eq. (28), the optical field interacts with 



all possible transition dipoles which couple the single- 



exciton states /i^^ 



the bi-exciton states 

x,xx 



and the 



a6' * 

single- to bi-exciton states /i^^^^ = M^^'^ • No permanent 
dipoles are present in this ensemble, i.e. /i^^ = = 0. 
In general, all these transitions are allowed due to the 
Coulomb scattering processes. 

Next, we employ the density matrix formalism to in- 
clude the dissipation processes due to the coupled phonon 
bath. Within this formalism, the dynamics of interest are 
fully described by the Liouville equation: 



(30) 



where the time-dependent density operator is a 2 x 2 
block matrix 



m = 



P^^'^{t) p^^ 



it) 
it) 



(31) 



containing single-exciton p^{t) and bi-exciton p^^(t) com- 
ponents, and coherences between single- and bi-exciton 
states p^^'^(t). The specific form of the relaxation term. 



T^p, in Eq. (30) depends on the specific exciton-phonon 



interaction model. 



A. Phonon-assisted relaxation model 

To describe the phonon-assisted dynamics, we assume 
that the phonon bath has a continuous spectral density, 
and there is no phonon bottleneck.' ^^W^TI explicit form 
of the spectral density depends on the environment model 
with adjustable parameters such as spectral widths and 
electron-phonon coupling strengths. The simplest model 
which can be employed in our case is the model of single- 



and bi-exciton states linearly coupled to the phonon co- 
ordinates {qa}a=i,2,3,...- The related Hamiltonian is 



Hp 



Hi + Hp, 



(32) 



where the exciton-phonon interaction term 

ab;a kl;a 

+ Y (l^«)^r/cfa + \xxk)Y^^.'^{xa\) , (33) 



contains the intraband single-exciton (bi-exciton) cou- 
pling matrix elements Y^^.^ {Yaj^a) ^"^^ phonon 
mode, and the interband coupling matrix elements Y^^.'^ . 
The connections between the former quantities and the 
e/ectron-phonon coupling constants from the many-body 
Hamiltonian are given in Appendix [B] The second term 
in Eq. ( 32 ) is the uncoupled phonon Hamiltonian whose 



form depends on the specific environment model. 

Assuming weak exciton-phonon coupling, we follow a 
standard projection operator method to eliminate the 
bath degrees of freedom resulting in the Markov approx- 
imation for TZp.^^ The basis set in which the equilib- 
rium (Gibbs) distribution, p, can be recovered as the zero 
eigenfunction of the relaxation operator, i.e. 7Zp = 0, 
is the quasipartide basis {|0 If =0,1,2,... formed by the 
eigenstates of the total projected Hamiltonian (Eqs. ([6|- 
64 Therefore, we consider the population relaxation 
dynamics in this preferred basis. 

After introducing the quasipartide energies, fkj^, and 
further using the interaction representation for the den- 
sity operator, i.e. p^/^/(t) = e~*^fc^£g(t) with 



we recast the Liouville Eq. (30) in the absence 



of the optical pulse {V{t) = 0) in the quasipartide basis. 
This results in the Redfield Equation: 



where IZ 



is the relaxation tensor.!^ 



The interaction representation allows us to apply the 
so-called secular approximation, eliminating the rapidly 



oscillating terms containing 



'a- 



-UJ 



the remaining Redfield tensor components,!^ 



^ 0. Respectively, 



R 



(35) 



correspond to uncoupled equations for the population re- 
laxation and coherence dephasing. Eqs. (35) and (36) 



contain the population relaxation and pure dephasing 
rates: 

% = ^E%V^f'«-;«'^-'KT)' (37) 
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(38) 



respectively. Here, C'^c^' (^) denotes the Fourier trans- 
form of the phonon correlation function, C^^' (r) = 
^^HpT-^^-iHpT-^^^^^^ and it has both real C'^^,{uj) and 

imaginary C^^, (u) parts. The explicit representation for 
the correlation function depends on the chosen form of 
Hp, i.e. on the specific relaxation model. 

The products of the off-diagonal quasiparticle-phonon 
coupling constants entering the population relaxation 
rate (Eq. ( [37[ )) can be expressed in terms of the exciton- 
phonon matrix elements entering Eq. ([33]) as: 



ll'rr' 



Finally, the matrix element of the diagonal quasiparticle- 
phonon coupling determining the pure dephasing rate 
(Eq. ^) is 



(40) 



In Eqs. (39) and (40), A^^(co'^) = lim^^o i'es{G^r(%)} 



is the transition amplitude given by the Green function 
residue in the limit of infinitesimal imaginary part, 7, of 
the poles. ''^ Here and below, we use the convention that 
the summation indices for single- and bi-exciton states 
(particularly those in Eqs. (39) and ([4Q|)) run over all 
single- and bi-exciton states, unless the superscripts x or 
XX constraining their range are used (e.g. in Eq. (33) for 

The dependence of the population relaxation rate 
(Eqs. ( [37[ ) and (|39|) on the transition amplitude indi- 
cates that these exciton scattering processes are involved 
in phonon-assisted cooling. Some of them, as we demon- 



strated in Sec. H C change the multiplicity of the initial 
states, and, therefore, can be considered as generalized 
impact ionization and Auger recombination processes. 
Accordingly, we argue that the phonon-assisted cooling 
and impact ionization dynamics giving rise to CM are 
generally coupled. However, we demonstrate in Sec. [IV| 
that the intraband phonon-assisted cooling and the in- 
terband impact ionization and Auger recombination pro- 
cesses can be decoupled in the limit of weak Coulomb 
coupling. 

By applying the secular approximation, we signifi- 
cantly simplify the description of the phonon-assisted 
dynamics. However, the validity of our approximation 
for NCs is based on the following delicate interplay be- 
tween the number of quantum states and their energy 
separations: In the region of high DOS, some closely ly- 
ing levels may have — uj^'^' ^ 0, potentially leading 
to the breakdown of the secular approximation. On the 
other hand, we expect that, due to the same high DOS, 



(a) 



(b) 



1^1X61 



(c) 



IC)(CI 
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-t' 






h 


G{ti) 






-ti) 


P 



Q{T-t') r 



gUu) 
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^pmi^ ) 













^pm{t'-ti) \g){g\ 



^pm(t'-ti) \9){g\ 



FIG. 5: Double-sided Feynman diagram representation of the 
nonequilibrium density operator, y5(r), prepared by the pump 
pulse and propagated during delay time, r. The times ti and 
t\ are the integration variables, and p is the equilibrium den- 
sity operator. p(r) can be partitioned into two components 
reflecting the contributions of (b) the quasiparticle coherences 
and (c) the quasiparticle populations. 



there are enough terms in the sum of Eq. ( 34 ) containing 
these slowly-oscillating phases to cancel out their con- 
tributions. Therefore, the chosen secular approximation 
must be validated using numerical simulations, for spe- 
cific materials. 



B. Photoexcited population dynamics and QE 

The central quantity describing CM is QE which can 
be calculated as 



QE 



2iV..(r)+7V.(r) 
^'..(r)+A^.(r) ' 



(41) 



where Nx{r) = trp^{r) and Nxx{'t) = trp^^{r) are the 
total non-equilibrium single-exciton and bi-exciton pop- 
ulations, respectively. Their dependence on the delay 
time, r, measured from the center of the pump pulse, 
allows one to calculate both the QE due to the photo- 
generation event and the total QE after the population 
relaxation. The latter is typically measured in optical 
experiments. 

The calculation of QE requires the solution of the Liou- 
ville equation given by Eq. ( 30 ) for single- and bi-exciton 
populations, p^. The solution of Eq. (30) can be ob- 



tained by using the fact that the coupling between the 
optical field and the exciton states is weak compared to 
the transition energies. This results in the second order 
perturbation expression associated with the double-sided 
Feynman diagram presented in Fig. [s] (a) The resulting 
population matrix element is 



lohh rorir2 
l~^liloPlorol~^rir2 

POO 

XI dt dti Qss,l2r2{^ ~ ^ 

) ^0 



(42) 



9 



where Gir(t) is th e matrix element of the projected prop- 
agator (Sec. |lIB ), and Gii^isni^ — t ) is the matrix ele- 
ment of the Redfield equation (Eq. ([34|) Green function 
transformed to the bare single- and bi-exciton basis. Fi- 
nally, pir is the matrix element of the equilibrium density 
operator 



\xo){xo\ 



(43) 



k>l 



with 



0,k 



\x 



\x 



E 

k>l 



TrX,XXTrXX,X 



(44) 



Equation (42) can be used for numerical calculations 



of the QE. This expression is quite general. Its form does 
not assume that the secular approximation (allowing the 
decoupling of the coherence and population relaxation 
dynamics) is used. For further analysis, we partition 
the contributions to Ps(r) induced by the optical exci- 
tation of the quasiparticle populations and coherences. 
For this purpose, we represent the time-dependent ma- 
trix elements of the projected propagator as 



Glrit) 



^ A;r(wj)e" 



(45) 



where the quasiparticle complex frequencies uj^ 
are the poles of Girioo)^ and Kiriuo^) = Tes{Gir{oJ^)} are 
the complex transition amplitudes given by the Green 
function residue Eq. (45) clarifies the physical mean- 



ing of the latter quantity showing that this is a proba- 
bility amplitude for the transition between / and r states 
in the single- and bi-exciton basis associated with the 
propagation of the quasiparticle state, 

Substitution of Eq. (45) into Eq. (42) and partition- 



ing the quasiparticle coherence and population dynamics 
(secular approximation) allows us to recast ps into a sum 
of the two terms. 



(46) 



Here, the first term corresponds to the double-sided dia- 
gram shown in Fig. [5](b). It describes the contributions 
of the quasiparticle coherences: 



(47) 



'^2 



■ ^pm ) 7 



where fissoi^i) is the projection of the transition dipole 
moment between the quasiparticle ground and f-th states 
onto single-/bi-exciton states, 

^ssoiO = Y^sh{oJ^)Phl^pl^so^ (48) 

containing the matrix elements, pi^so^ of the equilibrium 
density operator (Eqs. (44 )).''^ Note that PssoiO mixes 
the interband and the intraband dipole transitions en- 
tering the optical interaction term of the Hamiltonian 
(Eq. ([28|), and determine all possible photogeneration 
pathways. 



The pulse self-convolution function in Eq. (|47|) is 



prrij 
I 

t 



dti 



(49) 



This function is weighted by the coherence between the 
quasiparticle excited and ground states oj^.^q = uj^. —ojo — 
7^.0, and by the excited state coherences characterized by 
cD^. = uj^. — —7^.^., where the dephasing rates, 7^.^ 
and 7^.^^. , are determined by Eq. (38). 

The second term in Eq. p6|, represented by the 
double-sided diagram in Fig. p] (c), describes the con- 
tributions of the quasiparticle populations 

Ir 

M^so(OMrso(^)^(%0 

Here, Air{oJE) = lim^^o Azr(^£), and ei'^) ^CC eei'^) 



denotes the quasiparticle population relaxation compo- 
nent of the Green function associated with the Redfield 
Equation. This Green function can be found in the stan- 
dard way by using the eigenstates and eigenvalues of the 
Redfieled operator (Eq. ^).^^^'^^^ If the high DOS 
does not allow the diagonalization of the relaxation oper- 
ator, then 77,5(0) should be considered as the initial con- 
dition for the numerical solution of the Redfield equation 
(Eq. ( 34 j).!^ Finally, the pulse self-convolution function 
in Eq. (50) simplifies to the form: 



2 roo poo 

I{Cj^q - cupm) = ^ J dt J (51) 

X COS [(cj^o - ^pm) h] £pm{t )£pm{t - h), 

where we neglect the population relaxation processes dur- 
ing the interaction with the pulse. 



The representation given by Eqs. (46)-(51), provides 



a connection with the sum-over-eigenstates representa- 
tion shown in Fig. [s] (b) and (c). In this representation, 
there is an additional Liouville space pathway contribu- 
tion to the single- and bi-exciton populations associated 
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with the propagation of the ground state wave packet.'^ 
This term can, in principle, contribute to the CM dy- 
namics if Coulomb coupling between the vacuum and bi- 
exciton states is strong enough to make the lowest excited 
state energy comparable with the thermal energy, ksT^ 
i.e. V^^'O 2Eg — ksT. Since, the latter condition is 
not satisfied in NCs where typically Eg ^ yxx^o ^ ksT^ 
we do not consider this pathway. 

Finally, one can expect that the contribution of the 
quasiparticle coherences (Eq. ([47|)) to QE can become 
negligible compared to the quasiparticle populations 
(Eq. ([50|). This could happen, since spectral widths of 
ultrafast pump pulse can excite a significantly large num- 
ber of states (Fig. [Ij^a)), whose phases entering Eq. (47) 
through Eq. (|49| add destructively. This assumption can 
be checked for specific materials through numerical eval- 
uation of the related terms. 



IV. LIMIT OF WEAK COULOMB COUPLING 

In this section, we consider the Exciton Scattering 
Model, developed in Sees. [TI| and |III[ in the limiting case 
of weak Coulomb coupling. This limit is important for 
applications and assumes that the Coulomb matrix el- 
ements between single- and bi-exciton states are much 
smaller than the energy differences between these lev- 
els and/or much smaller than the level broadening, i.e. 
^ak^ <C h {\uj^ — ^k^l^Ja'k^)- demonstrate be- 

low, both carrier photogeneration and population relax- 
ation dynamics can be described using no higher than 
second-order processes in the Coulomb expansion. For 
this purpose, we use Eqs. (46)-(51 ), with the Green func- 



tion components calculated in this limit. 



A. Time-domain Green function 

To find the Green functions, we, first, represent the 
single-exciton (n = x) and bi-exciton (n = xx) free prop- 
agators (Eq. ([l6l)) as 



i5. 



kl 



(52) 



where the complex frequency, o)^ = cj^ — ^7^, contains 
the k-th frequency, cj^, from the projected Hamiltonian 
(Eq. ([7|), and the related dephasing rate, 7^. 

If the interband Coulomb interaction is weak, the 
CM dynamics becomes dominated by the Born inter- 
band scatting represented by the first vertex diagram in 
Fig. [3](c). According to Eq. (26), the scattering matrix 



elements in the Born approximation become 

TiX,XX / \ / • i-\—1t rX,XX 



(53) 



Furthermore, the leading contribution to the even-order 
scattering matrix (Fig.[3](a) and (b)) comes from the self- 
energy, and according to Eqs. (24) and (25), its single- 



and bi-exciton components become 



TrX,XXTrXX,X 
k ^k,b 



k>l 



TrXX,XTrX,XX 



o>0 



^2(0 



(54) 
(55) 



respectively. 

To calculate the time-dependent Green function, we 
substitute Eqs. (I52l)-(l55l) into Eqs._j2Tl)-(l23l). Further 



use of Fourier transformation (Eq. (llSfj leads to the fol- 
lowing expressions: 

G:,{t) = .5„fce-*^«* + Al, (e-*^«* - e-'-i^*) (56) 

k>l 

Glfit) = Skie-'^^^^* + Alf (e-*"^'* - e-*^?'*)(57) 

EkXX,X \X,XX -iu'^t 



a>0 



G 



X,XX f,\ kXjXX ( 

a,k y'^) - ^^a,k y 



XX j.\ 

' )■ 



(58) 



Here, the shorthand notations for the renormalized com- 
plex single- and bi-exciton quasiparticle frequencies 



(59) 
(60) 



are used, respectively. They contain the following self- 
energy corrections 



^rX^XX-rrXX^X 

^a,k ^kA 



E 

a>0 



■xrXX^X-xrX^XX 

k,a a,k 



(61) 
(62) 



Finally, the transition amplitudes in Eqs. (56)-(58) are 



\x 



A XX 
^^kl 



(1-ME 



k,b 



k> 



(1 - 40 E 



k,a a,l 



a>0 



(63) 



(64) 



Al 



-jXX 



(65) 



This representation for the time-domain Green func- 
tion (Eqs. (56)-(65)) is accurate up to second-order terms 
in the interband Coulomb interactions. In the following, 
the above expressions are employed to provide the lead- 
ing contributions to the single- and bi-exciton photogen- 
erated populations and to derive a set of rate equations 
for the population relaxation. 
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B. Single- and bi-exciton photogeneration 

The use of the Green functions represented by 



Eqs. ([56|)-([65| together with Eqs. (|46|)-([5l| results in 
the following form of the photo-generated single-exciton 
population: 



(66) 



where zeroth-, first-, and second-order terms describing 
the contributions due to the optically prepared quasipar- 
ticle populations are: 



oo \ ^ x,xx\xx,x X n-f-x 



k>l 



(67) 
(68) 

(69) 



k>l 



2Re {fioaKbf^bO + ^la^lhKo) 



b>0 



+ 2Re J2 f^oaKTt^lf^tr } ^(^a - ^P^l 

kl>l 

respectively. They contain the pulse self-convolution 
function (Eq. (51)) which is resonant only at single- 
exciton quasiparticle frequencies. 

In contrast, the quasiparticle coherences contributing 
to contain both single- and bi-exciton resonances 



(70) 



6>1 



X X{UJ^ — UJpm'-, ~ ^pm) 



XX 

k 



k>l 



6>1 



^pm ) 7 



entering the pulse self-convolution function (Eq. ([49|)). 
Although Eqs. (67)-(70) are important for numerical cal- 



culations of the QE, we do not discuss the scattering 
pathways associated with each term, since these path- 
ways carry no information about the CM dynamics. 

The photogenerated bi-exciton population in the weak 
Coulomb limit according to Eqs. (46)-(51 ) and Eqs. (56)- 
(|65| is 



XX 

Pk 



c(2) 



K2) 



^/c ^cr"~'^ (71) 

where the quasiparticle population contribution is 



xx{2) 



E 



(72) 



a>l 1>1 

and the quasiparticle coherence contribution is 



ci2) 



= -2Re PVKTKr^w (73) 

iay^b)>l 

- 2i?eEMSa*A::^ 



•^xx 
,k 



a>l 



EkXX,X'^X * I \ ^ ,XX^ kXX,X 
^^k,b PbO + Z^Pkl ^H,0 

b>l 

X -UJpm^UJ^ 



1>1 
- UJpm ) • 



Note that Eqs. (|71[)-(73) contain only second-order 
Coulomb terms. 

Equation ( [72| ) has a clear physical interpretation, il- 
lustrated in Fig. [6) where two interfering photogenera- 
tion pathways can be distinguished: The first pathway 
is shown in panels (a) and (b) and both the vacuum 
to single-exciton dipole transition (//Qa) interband 
Born scattering (A^'^^). The product, PoaK\T ^ 
latter quantities enters the first and second summations 
over single-exciton index a. These sums describe the re- 
distribution of the bare single-exciton oscillator strength 
between the quasiparticle single-exciton (panel (a)) and 
the bi-exciton (panel (b)) resonances. The second path- 
way, shown in panel (c), is represented by the last sum 
over the bi-exciton index, Here, the optical transition 
to the quasiparticle bi-exciton resonance is a combina- 
tion of the Born scattering event between the vacuum 
and a bi-exciton state (Aq'^^) and the bi-exciton intra- 
band transition (/i^f). Comparison of Eqs. (72) and (73) 



shows that the latter contains the interference of similar 
scattering pathways. 



C. Population relaxation 

To derive a set of rate equations for the popula- 
tion relaxation, we first represent the Redfield equation 
(Eq. (|34|) in the bare single- and bi-exciton basis. In 
this representation, populations and coherences are cou- 
pled. We eliminate the coherences and obtain a memory 
kernel that depends on the interband Coulomb coupling. 
We further apply the Markov approximation to the kernel 
based on the main assumption that the Coulomb inter- 
action is much smaller than the line width arising from 
the pure dephasing processes. This procedure results in 
the following set of rate equations: 



■X _ _ i,X,XX / X _ xx\ 

Pa — '^a,m \Pa Pm ) 



(74) 
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(a) 

Exciton Bi-exciton (b) 

Exciton Bi-exciton (c) 

Exciton Bi-exciton 

band band band band band band 

|x) \xx) \x) \xx) \x) \xx) 




FIG. 6: Bi-exciton photogeneration pathways in the weak Coulomb hmit. Panels (a) and (b) show the two components of the 
pathway involving the vacuum to single-exciton dipole transition (/j^) and the interband Born scattering (^A^'^"^). In panel (a), 
the intraband transition is in resonance with the optical pulse (hupm) but the final bi-exciton energy is distributed around 
hujpm according to the non-zero components of A^^^^ , Panel (b) describes the opposite situation, where the single-exciton is 
virtual and final bi-exciton state is in resonance with the optical pulse, (c) The pathway containing the production of virtual 
bi-exciton states due to the Born scattering from the exciton vacuum (A^'"^"^) followed by the intraband dipole transition /j^^ . 
Here, the final bi-exciton state is in resonance with the optical pulse. 



■XX 

Pk 



b 

-EK'ripr-pt) 

b 

EZ-pXX XX -pxx ^xx\ 



(75) 



Here, the interband scattering and the intraband 
phonon-induced population cooling are described by dif- 
ferent terms indicating that in the weak Coulomb limit 
these two processes are uncoupled (Fig. [t]). Specifically, 
the first term in the r.h.s. of each equation describes the 
interband population transfer due to both the impact 
ionization and the Auger recombination processes. The 
related population transfer rate, arising from the Markov 
kernel appearing in the coherence elimination, is 



the form of the population relaxation rates does not 
change. As a result, their expressions are 

Tab = ^ E ^a1,;a>rM' K " ^t), (77) 

r^r = ^E^'Sr«nT«'C'a.'('^r-c^r), (78) 



where the phonon correlation function, C^^'{uj)^ is de- 
fined in Sec. |III A} and the exciton-phonon coupling con- 
stants, Y^^.^ and Y^^^^^ are defined by Eqs. (B3) and (B4), 
respectively.^^ Similarly, the following expression for the 



pure dephasing rate immediately follows from Eqs. (38), 
(|40|), (l56|), and (Isrl) 



l,x,xx 



T7'X,XX 

^a,k 



jXX 



(76) 



where is the pure dephasing rate. 

The second terms in the r.h.s. of both Eq. (74) and 



(75), describes phonon-assisted cooling. The entering 



population decay rates can be obtained from the gen- 
eral expression given by Eq. (37) and (39), where we 
set Kbii^x) = ^ab, ^kii^xx) = and Aak{(^x,xx) = 
0. These transition amplitudes arise from the zero- 
order Coulomb terms of the Green functions given by 
Eqs. ([56])-(58). The remaining second-order transition 



amplitudes are dropped, since together with the exciton- 
phonon couplings their net contributions to the rates be- 
come negligibly small. 

Since, during the transformation from the quasiparti- 
cle representation back to the bare single- and bi-exciton 
states, one has to keep only zeroth-order Coulomb terms. 



ia,n 



/ J V aa\a. 



(79) 



(y-^ , -r^^ ,)c' ,(0). 



Finally, we outline the computation of the QE in the 
weak Coulomb limit. First, Eqs. ( 66 )-( 73) are evaluated 
to find the initial conditions, p^(Oj and p|^(0), for the 
population relaxation. These density matrix elements 
can also be used to obtain the QE (Eq. (41)) associ- 



ated with the photogeneration processes. Next, starting 
with the latter boundary conditions, the wavepackets, 
Pa{t) and p^^{t)^ should be numerically propagated to 
the bottom of the single- and bi-exciton bands accord- 
ing to Eqs. and (75) with the parameters defined by 
Eqs. ([76l)-([79|^ This provides input for the determina- 
tion of the total QE. 



13 



Exciton 
band 




(a) 



FIG. 7: Population relaxation in the weak Coulomb limit 
consists of uncoupled interband Auger recombination and im- 
pact ionization processes with rate k^^^^^ and the intraband 
phonon-induced cooling with rates and F^^. 



V. DISCUSSION 

In this section, we discuss the relation between our pro- 
posed Exciton Scattering Model and three earlier mod- 
els: the Coherent Su perpos ition Model, the Direct 
Photogeneration Model,^^^Ell ^^d the Impact Ionization 
Model^^^^. 



A. Coherent Superposition Model 

The Coherent Superposition Model is the limit of our 
Exciton Scattering Model, in which only the two states 
\xc) and \xxc) are coupled by the Coulomb matrix ele- 
ment V^'^^, and the two states \xu) and \xxu) are de- 
coupled as shown in Fig.[8](a). There is also no coupling 
to the vacuum state. The \xc) and \xxc) states are as- 
sumed to be almost degenerate, i.e. huj^^ ~ leading 
to the strong interaction condition h\uj^^ — cj^l <C V^^'^. 
The details of the calculations of the scattering matrix 
components, the single- and bi-exciton Green functions 
in the framework of the Coherent Superposition Model 
are given in Appendix[D] According to these calculations, 
strong Coulomb interaction corresponds to the splitting 
between coupled states, and formation of the quasipar- 
ticle states |±) with energies, huj± (Eq. (D5)), as illus- 
trated in Fig. [s] (b). Another assumption used in the 
calculations is that the splitting, cj^ = cj^ — uj- , signif- 
icantly exceeds the quasiparticle level broadening. These 
assumptions in our new notations reproduce the model 
proposed in Ref. i46] 

The authors of Ref. Hfi' used the phonon-assisted relax- 
ation model containing uncoupled intraband relaxation 
pathways for single- and bi-excitons resulting in inde- 
pendent cooling within each manifold. We argue that 
these relaxation pathways are coupled since the Coulomb 
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F..._ 
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FIG. 8: Level diagram for the Coherent Superposition model: 
(a) Bare single- and bi-exciton state representation where two 
states \xc) and \xxc) are coupled by Coulomb matrix element 
y^'^^, and two states \xu) and | 

xXu) are uncoupled. The 
transition dipole, /x, couples the vacuum state with the upper 
single-exciton states, (b) Quasiparticle representation. Scat- 
tering processes correspond to the optical transitions with 
/i+ = y/A^/j. and /i_ = y^A^/i. In a short-pulse limit, the 

pump spectral width exceeds huj-i , and both the coherence 

and populations of |it) states are included. Their projections 
back to the bare single- and bi-exciton populations are given 
by Eqs. (83) and (84). The p opul ation relaxation pathways 

-(86) are shown by the arrows 



(83) 



with rates entering Eqs. 
pointed down. 



assisted processes should be accounted for. Furthermore, 
we point out in Sec. |III A[ that the relaxation equations 
should reproduce the equilibrium distribution function 
not for the bare single- and bi-exciton states but for the 
quasiparticle states. Another assumption used in Ref. [46l 
is that the dephasing rate between coupled states is fully 
determined by the population relaxation processes. How- 
ever, the pure dephasing time in NCs is estimated to be 
several orders of magnitude shorter than the population 
relaxation time.^^ In this case, the dephasing rate is 
totally due to the pure dephasing, indicating separation 
of the timescales for the coherence and population dy- 
namics. 

The population relaxation rates from the quasiparti- 
cle states to the uncoupled single- and bi-exciton states, 
shown in Fig. [s] (b), immediately follow from Eqs. (37) 
and (39) with the t ransition amplitudes A^ and A^^ 
(Eqs. ([D9l) and ( |D1QD ) inserted. The rates are 



= ^E[^nc;j'A±C'.K-C.±), (80) 

a 

r-,± = pEKcTJ'A|-aK--^±), (81) 



interaction is strong. As a result the interband phonon 



where Y^^.^ and Y^^^^ are the exciton-phonon interac- 
tions connecting the uncoupled and coupled states. Note 
that these expressions reproduce the uncoupled single- 
and bi-exciton relaxation rates only if the Coulomb inter- 
action is weak, h\uj^^ —uj^\ > F^^'^.l^ Population trans- 
fer also exists between the quasiparticle states. However, 
for the sake of simplicity, we drop this pathway. The 
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derivation of the pure dephasing rate, 7+-, between the 
quasiparticle states is discussed in Appendix [P] 

The Coherent Superposition Model consists of only 
one optical transition induced by the dipole moment, 
/i, between the vacuum and coupled single-exciton state 
shown in panel (a) of Fig. [s] The multiple-scattering 
processes redistribute the oscillator strength so that both 
quasiparticle states become optically allowed as shown in 
panel (b) . To observe the oscillations of the pump-probe 
signal (bleach) predicted by the Coherent Superposition 
Model, the pulse duration should be less than the de- 
phasing time and spectral widths of the pulse should ex- 
ceed the level splitting, This condition is satisfied 
in the so-called impulsive limit in which the pulse self- 
convolution function becomes frequency-independent^^ 



= 



(82) 



-2-2 



Here, is the amplitude of the pump pulse, and = 
^fiiT^rn is the effective pulse duration. 

To find the time-dependent populations of the cou- 
pled states prepared by a short pulse, we use the ex- 
pressions for the corresponding density matrix compo- 
nents (Eqs. |D14l ), dmi] ), ( [D15l), and ( [D2T] )) obtained 
in Appendix |D| along with Eq. (82). This results in the 



following populations of coupled single- and bi-exciton 
states: 



+ 2A2A^A^cos(a;+_r)e-^+-^, 



(83) 



(84) 



l2 



- 2A^A^A^cos(u;+_r)e-^+-^, 

respectively. The populations of the uncoupled states 
as a function of delay time, r, can be obtained from 
Eqs. (|D16|) and (IDITI) together with Eq. (I82|): 



-(r,,c+r..,e)^ 



V -A? 

^ XX, 



■ xx,^ 



(85) 



(86) 



^/fi is a dimensionless parameter; A|, 



Here, A. — ipm/ v^j.j.j.xv^j.j.ijj.v^j.j.j.v^^ij j^ouj.t^j.j.xv^uv^x, ^^^7 

A|'^^, and A|^ are the transition amplitudes defined in 

Eqs. ([D9l)-([DrT|. 

The oscillations predicted by the Coherent Superposi- 
tion Model exist only between coupled states, and have a 
frequency equal to the quasiparticle level splitting, cj^ 



Specifically, they are described by the second term in 
Eqs. (83) and (84). Since there is no oscillating term in 
the populations of the uncoupled states, oscillations of 
the bleach can be observed in experiments in which the 
probe directly monitors the time evolution of and/or 
p^^. In fact, the expression for bleach, given by Eq. (7) 
in Ref. 46, contains non- vanishing contributions from the 
coupled bi-exciton population. 

Finally, we compare the QE of the photogeneration 
event, and the total QE after population cooling. The 



photogeneration QE can be determined from Eqs. (41), 



(83), and (84) on timescales longer than the dephasing 



occurs but shorter than the population relaxation. This 
quantity is: 



12(V^'^^)2 + {uj^^ -uj^y 

^Yx,xx^j2 ^ (^^xx -uj^Y 



(87) 



It contains no relaxation parameters. As the strength 
of the Coulomb coupling increases, QE^ approaches its 
maximum value of 3/2. 



The total QE can be determined from Eqs. (41), (85), 



and (86) in the limit in which time is longer than the 
typical population relaxation time. This results in the 
sum, 



2r 



Q^> = ^A|- 



xx,^ 



- x,^ 



where A| determines the probability of optical excitation 

for each quasiparticle state, ^ = ±, and the ratio of the 
population relaxation rates gives the maximum QE asso- 
ciated with each of the states. (Fig. [s] (b)). We empha- 
size, that this expression for QE accounts for the relax- 
ation pathways mixing the coupled and uncoupled states 
of different multiplicities. The latter contributions have 
not been considered before and are expected to be non- 
negligible. As the Coulomb interaction increases, QE^ 
approaches its maximum value given by the ratio of the 
total population relaxation rates (2r^^ +r^)/(r^^ + 



where F^ 



and F^ 



■ XX, -\- 



Regardless the differences in the relaxation models, sim- 
ilar maximum values for QEy were obtained in Ref. [46 
for increasing Coulomb interactions. 



Direct Photogeneration and Impact Ionization 
Models 



The Direct Photogeneration Model assumes weak 
Coulomb coupling between single- and bi-exciton states. 
Therefore, to find the bi-exciton generation rate, we 
should begin with the expressions for the photoin- 
duced bi-exciton population (Eq. ([tTI) 



-(73)) derived in 
Sec. |IVB[ An additional assumption of the model is that 
the pump pulse is much longer than the dephasing times, 
resulting in the so-called continuous wave (CW) limit. 
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The pulse self- convolution function for quasiparticle pop- 
ulations in this CW limit becomes proportional to the 
Lorentzian line-shape function: 



27r 



El 

a>l 



ttX^XX 



he: 



ET) , (93) 



i9 



(0) 



■pm ^pml^g 



t9 {uj, 



(89) 



where is the pump pulse electric field amplitude, and 
^pm = \/^^pm is the effective pulse duration. The pulse 
self-convolution function associated with the quasiparti- 
cle coherences vanishes in CW limit. 

Next, we introduce the bi-exciton generation rate as 
^xx = X]/c>i Pk^ /^pm where the limit of 7^0 ~^ should 



respectively. Comparison of these expressions with 
Eqs. (1) and (2) from Ref. 49, leads to the conclusion that 
the Impact Ionization Model is just the relaxation compo- 
nent of our model in the weak Coulomb coupling regime 
(Fig.[7|). Usually, the initial condition for the impact ion- 
ization dynamics is taken to be only the photogenerated 
5m^/e-exciton population (Eq. (67)). We argue that the 
bi-exciton population described by Eqs. (71 l-ffs]) should 



be taken. This, according to Eqs. (71) and (73), cor- 



responds to the following expression for the bi-exciton 
generation rate: 



jl 

k>l a>l 



XXX, X X 



5 {El - hujpm) 



k>l 



EA XX jX 
^^k,a 



PaO 



E,,XX A XX,X 
Pkl ^h,o 

1>1 



(90) 



also be included, since this contribution is of the same 
order of magnitude. 

The discussion above shows that the Direct Photogen- 
eration Model and the Impact Ionization Model comple- 
ment each other, and are a particular case of our more 
general Exciton Scattering Model in the weak Coulomb 
limit. Specifically, the Direct Photogeneration Model de- 
scribes the primary photoexcitation process involving a 
pump pulse that is longer than the dephasing time and 
shorter than the inverse relaxation rates of impact ion- 
ization (Eq. ([92])), Auger recombination (Eq. (|93|) and 
phonon-assisted decay (Eqs. ([77|)-([78|). The Impact Ion- 
ization Model describes photogenerated population re- 
laxation with the initial conditions given by Eqs. (90)- 



and E^^ 



huj^^ are the 



by Eq. (65), and E^ = huj[ 

single- and bi-exciton energies, respectively.^^ The lead- 
ing term in the single-exciton generation rate can be eas- 
ily obtained using the same approach: 



Here, theinterband transition amplitude, A^^ , is given ^ However, we argue that a systematic computational 



h 



^pm), (91) 



a>l 



and it coincides with that given in Ref. [27l 

Equation (90) should now be compared with 
Eqs. (1) and (3) describing the direct bi-exciton gener- 
ation rate via virtual single-exciton states and via cou- 
pling to the vacuum states derived in Refs. 27 and 47, 
respectively. This comparison shows that the two con- 
tributions to the second term in Eq. ([90| multiplied by 
5 {E^ — E^^) and shown in panels (b) and (c) of Fig. |6j 
reproduce Eqs. (1) and (3) from Refs. 27 and 47, respec- 
tively. In summary, the weak Coulomb coupling limit 
of our Exciton Scattering Model in the particular case of 
CW excitation not only recovers the previously developed 
Direct Photogeneration Model but also predicts an ad- 
ditional contribution given by the first term in Eq. (90) 
and illustrated in Fig. |6] (a) as well as the interference 
of the previously studied pathways (Fig. |6] (b) and (c)) 
which is clearly seen in the second term of Eq. ( [9Q| ). 

The central objective of the Impact Ionization Model is 
the calculation of impact ionization and Auger recombi- 
nation rates, which can be easily obtained from Eq. (76) 
by taking the limit of 7^'^^ 0, and further performing 
the summation over the final bi- and single-exciton states 
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approach should follow from the weak Coulomb limit 
computational scheme given in Sec. \P7\ since this ap- 
proach contains additional contributions not considered 
before and also accounts for finite-time pulse excitation. 



VI. CONCLUDING REMARKS 

Currently, numerical implementation of our Exciton 
Scattering Model is a challenging task, since the calcula- 
tions involve a large (~ 10^) number of bi-exciton states. 
This poses difficulties for the full scattering matrix calcu- 
lations (Eqs. ([24|)-([26|)) associated with large computer 
memory requirement for the matrix inversion. These dif- 
ficulties can be overcome by noticing that in NCs, the en- 
ergy difference between most of the coupled single- and 
bi-exciton states is larger than the interband Coulomb 
interaction, and only a small number of these states are 
in resonance (degenerate). As a result, the expressions 
obtained in the weak Coulomb limit (Sec. IV) should be 
used to evaluate the contributions from well-separated 
states, and only contributions from the degenerate states 
need to be included in the multiple-scattering formal- 
ism. If the level broadening for the degenerate states 
in NCs exceeds the Coulomb coupling then the full com- 
putational scheme developed in Sec. |IV| should be used, 
including the degenerate states. 

The comparison of CM processes in both NCs and in 
bulk semiconductors is important for understanding the 
role of quantum size effects on QE. Therefore, we em- 
phasize that our proposed formalism is valid for CM in 
bulk semiconductors. The transition is simple: One has 
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to replace all the summations over the single- and bi- 
exciton indices as well as over the quasiparticle states by 
summations over their quasi-momenta and spin degrees 
of freedom. All matrix elements entering the calcula- 
tions can be represented in the quasi-momentum basis 
set. This representation will automatically impose the 
quasi-momentum conservation restrictions. 

To summarize, we have proposed the Exciton Scatter- 
ing Model which treats the two main processes of CM, 
photogeneration and population relaxation, on the same 
footing. Our model is valid in the neighborhood of the 
AET where the contribution of the higher-multiplicity 
exciton states (tri-exciton, etc.) can be neglected. Our 
model includes relatively large Coulomb interactions 
leading to multiple interband scattering events. The only 
restriction on the Coulomb interaction strength is that it 
should not mix the higher multiplicity states. Based on 
our general formalism, expressions determining the QE 
in the limit of weak Coulomb interaction have been de- 
rived. This limit is extremely useful for numerical calcu- 
lations for specific materials. Since the AET is sensitive 
to the material-dependent selection rules imposed on the 
Coulomb matrix elements and transition dipole matrix 
elements, its determination can be done through direct 
numerical calculations. As we demonstrated, our Exci- 
ton Scattering Model recovers three previously proposed 
models as limiting cases. By including additional mech- 
anisms of CM, our model provides a unified approach to 
the study of CM in NCs and in the bulk limit. 
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Appendix A: Interacting multi-exciton Hamiltonian 

After introducing the electron and hole creation (an- 
nihilation) operators cj^ (cn) and d}^ {dm)^ respectively, 
the many-body Hamiltonian describing the valence and 
conduction band electronic states in semiconductors can 
be represented as a sum of three components^ 

?ie^=?io+Hl+H2. (Al) 

In this expression, the first term 



V ^mnlk ^mnkU^m^l ^n^k 



mnlk 



-yv 

mnlk 



eeee f f 
mnlk^m^n^l^k 



hhhhj]j]j J 
mnlk ^k^^^n^ 



mnlk 



conserving number of quasiparticles describes those non- 
interacting electrons and holes which are characterized by 
the Hartree-Fock energies and e^, respectively. Here 
and below (Eqs. ( |A2[ )-( [A5l )), the superscripts e and h 
denote electrons and holes, respectively, while the corre- 
sponding indices run over all Hartree-Fock states in the 
valence and conduction bands. In Eq. (A2) the Coulomb 
matrix element is 



Vr 



mnlk 



(A3) 

where V{\x — y\) is the Coulomb potential and i/jn is the 
electron (hole) Hartree-Fock wave funct ion. 
The second term in the Hamiltonian (Al), 
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ilk 



describes the processes of creation or annihilation of a 
single electron-hole pair (in the presence of another elec- 
tron or hole state) which are referred to as Auger recom- 
bination and impact ionizatio n, re spectively. Finally, the 
last term in the Hamiltonian (Al) 



1 



^^2 - 2 X] ^mnlk^mdnClCk 
mnlk 

+ 2 ^ ^mnlk^m^liAd\ (^^) 
mnlk 

characterizes those processes which involve the simulta- 
neous creation or annihilation of two electron- hole pairs. 

Neglecting all possible charged states, we consider only 
the space spanned by all possible multiple electron-hole 
pairs (multi-excitons) 5* = 0n>o5'^ where is the 
exciton vacuum (filled valence band, empty conduction 
band), and S'^ is the exciton space of multiplicity n with 
the complete basis set constructed from non-interacting 
electron-hole states 



|eX) = nLi4<|0), 



(A6) 



where the generalized indices are defined as a = 
{ai, . . . , ttn} and h = {61, . . . , 6^}. In this representation 
the many-body Hamiltonian term Ho (Eq. (A2)) maps 
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on itself. As a result, one can define the eigenstates 
\X^) of 7^0 forming a complete basis set in and, re- 
spectively, the eigenenergies, hw^. The latter eigenstates 
describe bound n-exciton states in whose eigenen- 
ergies include the binding energy due to the electron- 
electron, hole-hole and electron-hole Coulomb correla- 
tions. The introduced bound n-exciton states are related 
to the non-ineracting electron-hole basis through the uni- 
tary transformation 

defined by the matrix {/7^_^}. 
The eigenstate equations 

Y^iKerlHo - tuv^a\epK)UZ,p, = 0, (A8) 



In Eqs. (A8)-(A12) indices a and b denote the single- 



and bi-exciton eigenstates, respectively. The rest of the 
indices describe carriers states. 



Appendix B: The exciton-phonon interaction 
Hamiltonian 



The many-body Hamiltonian accounting for the linear 
electron-phonon interaction in the electron-hole represen- 
tation i^ 



mna 
mna 



and 



"^{heefhgehino - hj^''\ephqerhs)U^^pq^, = 0, (A9) 



where the phonon normal modes are q = 
{^1, ' ' 'Qa ' ' ' QNph}^ arid the coupling constants are 
the following matrix elements 



determining the transformation matrix elements for the 
single- and bi-exciton states, have the following form: 



+ K 
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{VmWaWn) ^ r, S = 6, h; 



(B2) 



{{^gq^es - ^eq^gs) [{^hp^fr - ^hr^fp) (4 + 4) (^^^) 



of the force operator Fq, averaged over the Hartree-Fock 
electron and hole wave functions. In this Hamiltonian, 
we dropped the term, fnn-a^ai which has no contri- 
bution to the processes under consideration. 

Mo dify ing the multi-exciton basis set given by 
Eq. (A7), the electron-phonon interaction Hamilto- 
nian (Bl) can be projected on single- and bi-exciton 
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- {SgqSes - ^eq^gs) {^hp^fr - ^hr^fp) ^^T) U'^lqrs = 0. 

respectively. Using the transformation matrix defined by 
the secular equations Eqs. (AlO) and (All) the Coulomb 
matrix element entering Eq. (jsl) becomes 
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states resulting in Eq. ( |33[ ) where the intraband exciton- 
phonon coupling matrix elements are 
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The interband exciton-phonon matrix element also en- 
tering Eq. (33) is 
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/ J / J ^ a;gf^ m;pqrsJ mn;a 
gpr=e fqs=h 

X {^fs^grfqp-a ~ ^fs^gpfqr;a 
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Here, the transformation matrices and U^^ can be 



calculated according to Eqs. (AlO) and (All), respec- 
tively. 
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Appendix C: The exciton-optical field interaction 
Hamiltonian 



where we denote the eigenenergies (quasiparticle ener- 
gies) as 



The many-body Hamiltonian in the Hartree-Fock or- 
bital representation describing the interaction with the 
time-dependent optical field E{t) is!^ 

HoAt) = -mY.P^'nClc„ (CI) 

mn 

+ E{t)Y,P:tdldn 

mn 

~ ^(^) Pmn^rnC-n 
nm 

mn 

where the transition dipole moments are matrix elements 

= {rj^-d\rn), r,s = e,h; (C2) 

of the dipole moment operator d projected onto the 
field polarization direction, e, and further averaged 
over the Hartree-Fock electron and hole wave functions. 
This Hamiltonian has exactly the same structure as the 
electron-phonon coupling Hamiltonian (Eq. (Bl)) dis- 
cussed in Appendix [B| Therefore, in direct analogy with 
Eq. (33), one can immediately recast the former Hamil- 



tonian to the single and bi-exciton state representation 
given by Eq. (28). The intraband single-exciton (bi- 
exciton) transition dipoles /i^^ (^mn) entering Eq. ( [28| ) 
can be determined by replacing Y^^.^ (^rfm-a) l.h.s. 
of Eq. (|B3| (Eq. ([bI|) by /i^^ (^ni)- Determination of 
the interband transition dipoles requires the replacement 
in the l.h.s. of Eq. (|B5|) hyjj^^. Also ah f^^.^ 



yx,xx 



(r, 5 = e^h) in the r.h.s. ot Eqs. (B3)-(B5) should be 
replaced by P^^ (r, s = e^h). 



Appendix D: Coherent Superposition Model 

Coherent Superposition Model is the limit of the Exci- 
ton Scattering Model including only two states \xc) and 
\xxc) coupled by Coulomb matrix element V^'^^ and two 
uncoupled states states \xu) and \xxu) characterized by 
the following Hamiltonian: 

^+ ^ {\Xa)huJ^{Xa\^\xXa)huJ^''{xXa\) (Dl) 
a=c,u 

+ \Xc)V'^''^'^{xXc\ + \xXc)V'^'^''^{Xc\. 



According to Eqs. ([24]) and (26), the scattering matrix 
elements for the coupled states are 
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Following the procedure in Sec. |H C[ one finds that 
the corresponding time-dependent Green function com- 
ponents for coupled states are 



A^e-^'^+^ + A^e 



Gl{t) 
GT{t) 

G^'^^(t) = A^'^^e-*^+^ + A! 
which depend on the transition amplitudes 



e 

x^xx 



Ax,xx 
± 



= ± 



± 



± 



{uj± - 




(^+ - 


a;_)' 


{uj± - 




(cj+ - 


uj_y 


yx 


,xx 


h (cj+ 





(D6) 
(D7) 
(D8) 



(D9) 
(DIO) 
(Dll) 



Note that these quantities are real, since the Coulomb 
coupling significantly exceeds the level broadening. 

To find the single- and bi-exciton populations associ- 
ated with the excitation of the quasiparticle populations, 
we use the following relaxation equations for populations 
in the quasiparticle representation (Fig. [s] (b)): 
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The non-vanishing (for r > 0) components of the Liou- 
ville space Green function (Eqs. ( D12[ )) are: 
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(D13) 



■ x,± 



■ XX, ^ 



(l-a±;±(r)), 
(l-a±;±(r)). 



where, the population relaxation rates Tx,± and Txx,± 
are given by Eqs. (80) and (81), respectively. For the 



sake of simplicity, we dropped the population transfer 
rates between the quasiparticle states. 

The time-dependent single and bi-exciton populations 
due to the excitation of quasiparticle populations ar e ca l- 
culated according Eqs. (48) and (50) where Eqs. (D9)- 



(Dll) and (D13) are substituted. For the coupled states 



these populations are 



^pm J •) 



), (D14) 



UpfYi 



), (D15) 
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and for the uncoupled states 

<{r) = /i^^a:,^-(r)A|X(c^f-c.,^), (D16) 

^=± 

<^{r) = /i2^^^|(r)A|X(a;^--c.,^). (D17) 

^"=± 

Here, the pulse self-convolution function is defined 
by E g. ( 51), and the population Green functions by 
Eqs. (|D13|). 



The Liouville equation for the coupled state coherence 
in the quasiparticle basis set is: 



(Dl^ 



where Y^^ , Y^^^^ , Y^J^^^ , Y^^ are the components of the 
exciton states coupled to the phonon mode, a. 



According to Eq. (47), the contribution of the quasi- 



particle coherences to the coupled single- and bi-exciton 
populations are: 



(D20) 
(D21) 



where u j^ = ujj^ — uj- — contains uj± determined 

by Eq. (D5) and the pure dephasing rate 7+_ which can 
be explicitly found by using Eq. (38) and the following 
quasiparticle-phonon couplings: 



where the pulse self-convolution function is given by 
Eq. (49). By taking into account that a^^'^A^^'^ — 



-AiA 



, one finds that c^(r) = — c^^(r). Obviously, 



l;^^A^^ (D19) cS(r) = (r) = 0. 
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The pump envelope and spatial phases do not contribute 
to the population dynamics and therefore are dropped. 
Using this representation for the transition amplitude, 
Eqs. ( 40 ) and ( 39 ) can be derived by straight forward trans- 



formation between the exciton and quasiparticle basis sets 
using the following relationship — ^ir{(^^) dis- 

cussed in Ref. 57 



We 



Eqs. (43)-(45) are obtained in the second order perturba- 
tion theory where the small parameter is V^'^^ /Eg 
also used the fact that in semiconductors NCs Eg ^ ksT. 
The poles are assumed to be of the first order which is 
the general situation for arbitrary Coulomb coupling. Two 
poles can coincide if we use the second order expansion of 
Eq. (21 ). This results in the renormalization of the energy 



in the exponential but the form of Eq. (45 ) will still be the 
same. 



Deriving Eqs. (47) and (48), we used the following relation- 



ship between the components of the equilibrium density 
matrix and trans itio n amplitude: Kir{ujQ) — pw- 
Note that in Eq. ( [50| for ns(0), the population Green func- 
tion becomes ^^^^(0) = J^^, and associated prefactor sim- 
plifies to [A5s(cJ<^)^^,^(0)Air(cJ^-)] = Ssl^sr- 

Formally, the expression for the interband population 

^ l^cxcx' ^ak-c, ^^a-oc' o^^' 



transfer rate F^L^^ 



xx\ 



can be introduced. However, in this expression Y/^^'^, ^ 

ka\a. 

only for the states which have h{uj^^ — ^a) > Eg. Since 
Eg significantly exceeds the phonon bath spectral widths. 

This follows from the observation that for y^'^^ 0, 
A^ = A^^ = 1/2 and A^ = A^" = 0. 

It is easy to show that for small Tpm the Gaus- 
sian pulse envelope can be approximated as Spm{t) — 

c(0) -^2/2x2 



27TrSpmS{t) and further substituted into 
Eqs. (|49| and (l51| to get the pulse self-convolution func- 
tions. 

The equalities require that Ca^uo) is a constant for u vary- 
ing within the range defined by y^'^^. 
This expression can be rigorously derived by using 
the Gaussian form for the envelope function 8pm{t) = 



(0) t2/2T 



Another way to derive Eq. ( 89 ) for a generic 



pulse envelope is to assume that Spm =const. during the co- 
herence integral over ti, and that the pulse auto-correlation 

function dt Spm{t )8pm(t - tl) ^ Tpm ^pm . 

The self-energy contributions to the energy are dropped 
because they correspond to fourth-order corrections. 
Notice that this equation does not account for the exciton 
filling factor effects since no more than a single- or bi- 
exicton states are excited per NC. In fact, the filling factors 
for the carriers forming the excitons are already accounted 
for by the exciton definition and enter into the Redfield 
equation implicitly. 



